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According to some recent analysis (M. Paczuski and M. Baiesi, Phys. Rev. E 69, 066106, 2004 
[1]) of earthquake data, aftershock epicenters can be considered to represent the nodes of a network 
where the linking scheme depends on several factors. In the present paper a model network of 
earthquake aftershock epicenters is proposed based on this scheme and studied on fractals of different 
dimensions. The various statistical features of this network, like degree, link length, frequency and 
correlation distributions are evaluated and compared to the observed data. The results are also 
found to be independent of the fractal geometry. 
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An earthquake is a complicated spatio-temporal phe- 
nomenon [2] which exhibits complex correlation in space, 
time as well as magnitude. Physically, earthquakes occur 
when the convective motion in the mantle cause sudden 
rupture and deformation of certain parts of the earth's 
crust that result in the emanation of energy in the form 
of seismic waves. In other words, earthquakes result 
from the interaction between stress concentration and 
fluid flow and have recently been the subject of consider- 
able interest in various fields including statistical physics 
because of its unique statistical features. 

Depending on the relative magnitude and position in 
the space-time sequence we can primarily classify earth- 
quakes into three categories of events [1-3] : 

(i) Intermediate or small amplitude precursor events that 
precede a main event, known as fore-shocks. 

(ii) Events of relatively large magnitude and impact, 
known as main shocks. 

(iii) Nearby smaller correlated events that follow a large 
seismic event or main shock, known as aftershocks. 

Earthquakes are best understood by studying large 
space-time correlations of many events instead of observ- 
ing isolated individual events and it has been suggested 
that earthquake is a critical phenomenon [4]. The strik- 
ing statistical features followed by earthquakes in general 
are the following: 

(a) The distribution of earthquake magnitude (m), is 
given by the Gutenberg- Richter (GR) law [5]: 

P(m) - 10-''™, (6- 1) (1) 

where P{m) is the number of earthquakes of magnitude 
TO in a seismic region. 

(b) The short time temporal correlation between earth- 
quakes follows the Omori law [6], which states that fol- 
lowing a main event, the frequency of a sequence of af- 
tershocks decays with time t as: 

iV(i)^t-",(«^l) (2) 



(c) The spatial distribution of earthquake epicenters form 
a fractal set with fractal dimension df [7]. 
Several analyses of real earthquake data reveal that 
earthquakes are a result of a dynamical many-body sys- 
tem that reaches a stationary critical state characterised 
by spatial and temporal correlations that follow power 
laws (eq.s 1,2) without any intrinsic time or length scales 
and hence they are related to self organised critical phe- 
nomena [8,9]. 

The Earthquake Network: 

In [1] the earthquake aftershock phenomenon has been 
visualised as a network. Networks are complex web like 
structures comprising of nodes, connected by links and 
such structures describe a wide variety of natural systems 
and have been studied recently with growing interest [10] . 
In [1] the nodes are the earthquake/aftershock epicenters 
and they are linked with a weightage given by the correla- 
tion between them. For the generation of an earthquake 
network, a metric is required which would quantify the 
correlation between two events or whether one event can 
be considered as an aftershock of another. This metric 
should incorporate the self similar statistical properties 
that unify earthquakes in general and is given by [1] : 



n^j = Ctri AmlQ-^"'' (3) 

where i and j represent events, i.e, earthquakes and/or 
aftershocks time ordered with i preceding j. 
The various quantities appearing in eq.(3) are as follows: 
iTii — Magnitude of the z*'* event; 

I — lij = Spatial distance between two earthquake epi- 
centers corresponding to the i*'' and j"* events; 
t = tij = Tj — Ti~ Time interval between the two events, 
with Tj > Ti, 

C = A constant depending on the overall seismicity in 
the region under consideration and 

riij = the expected number of events with magnitude 



within Am of rrii, occurring within the space-time do- 
main bounded by events i and j. 

Of all the earthquakes preceding j, there must be 

some event i = i* which is most unlikely to occur and 
for which iiij is minimum. Nevertheless since i* actually 
took place relative to j, inspite of being the most un- 
likely, so i* is the event most correlated to j. Hence the 
degree of correlation between any two earthquakes i and 
j is inversely proportional to riij and two events i* and j 
are linked if is minimum for i = i* . Thus a network 
is generated which is directed in time where, the nodes 
are the earthquake/aftershock epicenters, characterised 
by internal parameters which are magnitude, location 
and time of occurrence, and any two events are linked 
with a weight riij. In this paper our aim is to find out 
whether it is possible to construct a model which re- 
produces results comparable to the observed behaviour 
of the statistical features of some real earthquake data. 
We have also studied the dependence of the results on 
Euclidean spaces of different fractal dimensions. 

Observed results for real earthquakes: 
In [1], the catalogue maintained by the Southern Cali- 
fornia Earthquake Data Center (SCEDC) has been anal- 
ysed. The values of the various constants in eq.(3) that 
are used are b ~ 0.95, C = 10~^. 

The value of the fractal dimension used is df ~ 1.6 fol- 
lowing the result obtained in [11]. 
The analysis led to the following observations: 

(a) The distribution of the time intervals of the oc- 
currence of aftershocks, iy{t), shows a power law decay 
(vit) ^ ^ for t < tcutoff, c and K are constant in 
time and tcutoff depends on the magnitude m of the 
earthquake) which is consistent with the Omori law [6]. 

(b) The resulting network of earthquakes is scale-free, 
i.e., the out-degree distribution P{k) follows a power 
law decay, and for the SCEDC data it is found that 
P{k) ~ )fc-T with 7 = 2.0. 

(c) The link weight distribution N{nij) exhibits a power 
law decay with slope ~ 1. 

(c) The link length is defined in [1] as the distance be- 
tween the epicenter of an aftershock and its linked pre- 
decessor and its distribution L{1) shows an approximate 
power law decay. 

Modelling and Results: 
In our present work we simulate a network of earthquake 
aftershocks on Euclidean spaces with different values of 
the fractal dimension df. These are: 

(i) Continuous two dimensional Euclidean space {df = 
2.0) 

(ii) Percolation cluster on two dimensional lattice 
(df = 1.89) 

(iii) Backbone of a two dimensional percolation cluster 



{df = 1.6) 

(iv) Elastic backbone of a two dimensional percolation 
cluster (d/ = 1.10) 

After the generation of the network, we have studied the 
various statistical features, viz., time or frequency distri- 
bution, degree distribution, correlation distribution and 
link-length distribution. We compare the results for the 
four different types of Euclidean spaces with the observa- 
tions of real earthquake data and also examine whether 
the results at all depend upon the fractal dimensionality 
of the system. 

We have simulated the earthquake network on all four 

kinds of Euclidean spaces with different system sizes 
pertaining to technical limitations. An averaging over a 
maximum of 1000 configurations have been used in all 
cases. The network generation procedure is as follows: 

(i) Selection of nodes: For the continuous two dimen- 
sional Euclidean lattice, we take a unit square space 
where nodes occur randomly and each node is assigned 
coordinates Xi,X2 where < Xi < 1. 

In order to simulate networks on percolation clusters or 
its backbone or elastic backbone, we have taken a two 
dimensional square lattice where the occupation proba- 
bility is equal to the percolation threshold, pc = 0.592746 
[12]. Once it is checked that a percolation cluster (with 
df = 1.89) does occur, nodes are chosen randomly on 
this cluster. The same has been done for the backbone 
with df = 1.6 and elastic backbone with df = 1.10 after 
identifying these from the percolation cluster. 

(ii) Assignment of strengths: The selected sites represent 
the earthquake/aftershock epicenters or events and the 
epicenters are assigned strengths of magnitude to, ran- 
domly on a scale of 1 to 10. 

(iii) Calculation of Distance (kj): When the epicenters 
with their strengths are chosen, we calculate the Eu- 
clidean distance kj between any two events/epicenters i 
and j where i and j also denote their time of occurrence. 

(iv) Establishment of links: We start with the second 
event, j = 2 which is obviously linked with the only 
available preceding event, i.e, i = 1. For j > 3, events 
i = l,2,3,....,j — 1 are considered for the calculation 
of the correlation riij according to eq.(3) and node i is 
connected to node j for minimum . 

The values of the constants used are b = 1.0 and 
C = 10~^. While the maximum number of points used 
in the continuous 2d Euclidean space is n = 5000, for the 
percolation clusters we have used lattices of size L x L 
where L varies from 100 to 800. On the percolation 
clusters we have chosen from 500 to a maximum of 2000 
sites as epicenters. 

Statistical properties: 
After the links have been established, we have evaluated 
the different distributions defined earlier. The results for 
these distributions do not show any significant finite size 
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effect in any of the fractals considered. 

(I) Distribution of time interval between aftershocks: 
First we observe the nature of the distribution of time 
interval between events. This should follow the Omori 
law [6] , according to which the distribution should have 
a power law decay with an exponent ~ 1.0. We in fact 
obtain a power law fall with the exponent close to 1.0 for 
all the four fractals (Fig. 1). For the elastic backbone a 
cutoff region is observed beyond t ~ 100. 



(II) The out-degree distribution: The earthquake net- 
work is a directed network [1,10] where the links are 
directed in time. The in-degree, i.e, the number of in- 
coming links for any node is one, while the out-degree k 
of a node gives the number of aftershocks generated by 
that node. We observe that the number of nodes with 
out-degree k follows a power law decrease, P{k) ^ k^'^ 
with 7 ^ 2, for all four types of Euclidean spaces in 
agreement with the observations of [1]. Although the 
curves (Fig. 2) show occassional kinks, but an approxi- 
mate straight line can be fitted for all the four cases. 
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FIG. 1. The distribution of time interval between events. 
The distribution follows a power law decay in all cases with 
exponent close to 1.0 which is in agreement with the Omori 
law. A cutoff is observed beyond t > 100 for the elastic back- 
bone. 
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FIG. 2. The out degree distribution for earthquake lattice 
simulated on the four kinds of spaces lattice as mentioned in 
the figures above. In each case it closely follows a power law 
decrease with average slope around 2.0. 

(Ill) Correlation Distribution: The plot (Fig. 3) of the 
correlation distribution Mijiij ) also has a power law de- 
cay, with a slope ~ 2.6 ± 0.2 which deviates considerably 
from the observed value in [1] (this issue has been dis- 
cussed later). A cutoff region exists here for all the cases. 



tribution C{lij) for the network grown on the different 
spaces does not follow a power law decay except for the 
elastic backbone of the two dimensional percolation clus- 
ter, with slope ~ 1.2. For the other three cases, i.e., 
the continuous Euclidean lattice, the two dimensional 
percolation cluster and the backbone of the percolation 
cluster, the distributions show a power law decay with 
an exponential cut off given by, aa;"'' exp(— Aa;) where 
a, p and A are constants. The values of 0.1) and 
A ('^ 0.1) for the three cases are comparable; exact values 
are indicated in Fig. 4. 
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FIG. 3. The plots of the correlation distribution for the 
four types of fractals. Here the plot follows a power law de- 
cay in all cases with slope ~ 2.6 ± 0.2. 

(IV) Link length distribution: The link-length dis- 
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FIG. 4. The link length distributions for the four types of 
fractals. 
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To summarise the results, we have generated earth- 
quake aftershock networks on EucHdcan spaces with 
different fractal dimensions. The linking scheme follows 
the prescription of [1]. We have compared the distri- 
bution fimctions corresponding to the time intervals, 
out-degree, correlation and link-lengths from these net- 
works. The results partially agree with the real data, 
e.g., the time interval distribution shows a power law 
decay as given by the Omori law, the out-degree dis- 
tribution has an approximate power law decrease - in 
both cases the exponents also agree with the observed 
values quite well. The link-length distribution in [1] is 
approximately a power law while here it seems to have 
an exponential cutoff. For the correlation distributions 
the exponent value deviates appreciably from the ob- 
served one. While u{t), P{k) and £(Z) can be estimated 
independently without a network formalism, riij is a 
measure directly linked to the networking scheme. That 
our result for Af{nij) does not agree quantitatively with 
that of the SCEDC data [1] may be due to the fact that 
epicenters have been assumed to occur randomly while 
in reality there is expected to be a correlation. However, 
the exact values of Af(nij) (which depend on the pa- 
rameters chosen) may not be of much significance since 
the task is to link the epicenters with the minimum mj . 
The fact that we get good agreement for v{t) and P{k) 
and fairly good agreement for C{1) refiects this. On the 
theoretical side our studies show an interesting result 
that the behaviour of these distributions are almost in- 
dependent of the underlying spatial structure. We have 
varied the fractal dimensionality from 2 to 1.1 and we do 
not notice any appreciable change in the behaviour of the 
different quantities, either qualitatively or quantitatively. 
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